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Abstract. Within the framework of the scalar-tensor models of gravitation and 
by relying on analytical and numerical techniques, we establish the existence 
of a class of spherically symmetric spacetimes containing a naked singularity. 
Our result relies on and extends a work by Christodoulou on the existence of 
naked singularities for the Einstein-scalar field equations. We establish that 
a key parameter in Christodoulou's construction couples to the Brans-Dicke 
field and becomes a dynamical variable, which enlarges and modifies the phase 
space of solutions. We recover analytically many properties first identified by 
Christodoulou, in particular the loss of regularity (especially at the center), and 
then investigate numerically the properties of these spacetimes. 



1. Introduction 

The issue of the vahdity of the (weak version of the) cosmic censorship conjecture 
remains one of the most important open problems in classical general relativity. 
Roughly speaking, it says that physically admissible solutions to the Einstein equations 
should not contain naked singularites, that is, all singularities formed in physically 
reasonable scenarios of gravitational collapse should be surrounded by event horizons 
and, hence, cannot send signals to far observers at future null infinity. A precise 
formulation of the conjecture can be found in [28, 9], together with the properties 
required on a solution to qualify as a physically "reasonable" process of singularity 
formation. These properties concern the smoothness and genericity of the initial 
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conditions, and demand that the matter model undergoing collapse cannot form 
singularities of non-gravitational origin. 

Even though the conjecture is still far from being proven in general, no definitive 
counter-example has been found so far, neither in numerical simulations nor in 
analytical investigations. An important step forward in this respect was the numerical 
analysis of the threshold for black hole formation. After the pioneering work of 
Choptuik [10] it has become clear that it is possible to form a naked singularity by fine- 
tuning (smooth) initial conditions toward the vicinity of the threshold of formation of 
arbitrarily small black holes. It turns out that the process is dynamically controlled by 
an unstable exact solution -referred to as a critical solution- which, itself, does contain 
a naked singularity. The fine-tuning is required to compensate for the instability of 
this solution and achieve a continuous approach to that solution. In that set-up, there 
is no dynamical formation of naked singularities and, therefore, this analysis does not 
provide a genuine counter-example to cosmic censorship. (See [22] for a review.) 

In parallel to Choptuik's work on the critical collapse of a real massless scalar field 
in spherical symmetry, Christodoulou [7] studied the Einstein-scalar field equations 
from a fully analytical point of view. In a truly remarkable series of papers about the 
global dynamics of solutions to this system, he constructed a family of exact solutions 
parametrized by some reals fc,oi and showed that these spacetimes do contain a 
naked singularity in certain range of ai, provided < fc^ < 1/3. Later, in [8], he 
also established that these naked singularities are unstable under small perturbation. 
Christodoulou's work provided the first complete mathematical proof of the formation 
of a naked singularity under gravitational collapse. 

By construction, Christoudoulou's spacetimes are homothetic, that is, 
continuously self-similar and, therefore, do not contain any privileged scale (such as 
a horizon of finite size), and so cannot contain a (finite) black hole. Consequently, 
these spacetimes are a priori good candidates to contain naked singularities with a 
central point of infinite curvature, denoted by O; see [19]. The critical solution found 
by Choptuik also possesses self-similarity, but of a different type, known as discrete 
self- similarity. Since the symmetries are different, Christodoulou's solutions cannot 
"relax" to the critical solution, and actually the relation between the two solutions is 
unclear a problem that would deserve further study. 

The parameter k, whose origin is in the massless scalar field (which only enters via 
its derivative in the Einstein equations), is essential in Christodoulou's construction, 
as well as the key restriction 



Depending on the second parameter oi, an apparent horizon (rather than a naked 
singularity) is also possible in this range of k. For 



the future light-cone of O collapses to a line, which provides an example of a null 
singularity not preceded by an event horizon. For fc^ > 1 the spacetimes are rather 
pathological (see [4] and the Carter-Penrose diagram in Figure 4 of [15]). 

In all cases the parameter k generates a mild loss of regularity at the center, which 
makes the curvature to be continuous but non-differentiable before the singularity at 
O is formed. This implies that the past light cone of O is non-regular, hence the initial 
conditions are not completely regular. 

Note that, on the contrary, Choptuik's critical solution is smooth everywhere 
except at the central singularity at O and, as a matter of fact, the sole requirements 



< fc2 < 1/3. 



(1.1) 



1/3 < fc2 < 1 



(1.2) 
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of regularity and discrete self-similarity select a unique solution, at least locally in the 
phase space; sec [13, 14]. If the masslcss scalar field is taken to be complex then it is 
possible to construct a continuously self-similar solution which shares the regularity 
properties of the Choptuik spacetime, and also contains a naked singularity, but is not 
critical [17]. 

Our main objective in the present work is to investigate whether the relevance of 
k and the role played by the limiting conditions on k for the formation of naked 
singularities are "structurally stable" , that is, whether spacetimcs with the same 
features can be constructed with extended models in which families of solutions 
with an equivalent parameter are present. Indeed, the model we consider contains 
Christodoulou's model as a special case. Specifically, we work here within the scalar- 
tensor theory of gravitation and, especially, within the so-called Brans-Dicke theory. 

The model under consideration here effectively adds an additional scalar field to 
Christodoulou's system of equations, and makes k a dynamical variable, denoted by 
K. Interestingly enough, our analysis leads to the same range in order to avoid the 
pathological behavior referred to above, namely 

< < 1, (1.3) 

where now is the value of the field K at the past light cone of the singularity. 

We also show that, starting from an arbitrary initial value for the system under 
consideration dynamically evolves toward values K below the threshold 1. Therefore, 
the introduction of extra degrees of freedom allows to avoid the pathological spacetimes 
arising with the Einstein-scalar system. 

Note that Liebling and Choptuik [21] have numerically shown the presence of 
critical phenomena in the Brans-Dicke system, the critical solution being discretely 
or continuously self-similar (depending a coupling constant). Again, being completely 
smooth, such a critical solution is not related to the solutions we construct in this 
article. For further results, see [5,18,20,26]. 

The system under study is significantly more involved than the one studied 
analytically by Christodoulou [7] and, although we do follow and generalize several 
important steps in the construction therein, we eventually must resort to numerical 
investigations to reach our final conclusions. In fact, by relying on numerics, we arrive 
at a better understanding of the class of solutions and are able to construct explicit 
examples. 

An outline of this paper is as follows. In Section 2, we introduce the model of 
self-gravitating matter of interest, and we determine the general evolution equations 
under the assumption of radial symmetry. In Section 3, we impose the self-similar 
assumption and show that general solutions are parameterized by four functions of a 
single variable, denoted by x, which obey a system of ordinary differential equations 
(ODE). We construct solutions that are piecewise regular, with each piece separated 
by singular points across which careful matching is required. Specifically, in Sections 4 
and 5, we successively construct the interior and exterior part of the past light-cone 
of the singularity. Finally, in Sections 6 and 7 we describe our numerical strategy and 
present various results and conclusions. 
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2. Scalar-tensor theories 

2.1. Scalar-tensor gravity with scalar field 

Scalar-tensor theories of gravity are alternative theories of gravity which are physically 
strongly motivated and have a long history in the literature. The fundamental 
assumption of these theories is that the gravitational field is mediated by one (or 
more) scalar field(s) in addition to the standard tensor field (5^1,) of Einstein's general 
relativity. These theories satisfy the equivalence principle (since they are metric- 
based theories), but do not satisfy the strong version of the equivalence principle. 
The first theory of this kind was developed by Jordan [16], Fierz [12], and Brans and 
Dicke [6], and contains an additional parameter defining the coupling of the scalar 
field to the matter model. Later on, Bergmann [2], Nordtvedt [23], and Wagoner [27] 
extended this approach to a coupling via a function of the scalar field. Next, Damour 
and Esposito-Farese [11] introduced a generalization based on an arbitrary number of 
scalar fields. More recently, cosmological models based on the so-called f{R) gravity 
theories have attracted a lot of attention, which found applications in the study of 
relativistic stars [1] . These theories form a subclass of the scalar-tensor theories, and 
it is interesting to look for a better understanding of the corresponding spacetimes 
and, in particular as we do in the present work, to study the possible existence of 
spacetimes containing naked singularities. 

Specifically, we are going to investigate a generalization of Christodoulou's model 
when a a scalar field in coupled to a scalar-tensor theory of gravity for which the 
action reads (see [11, 25] for details): 

S = Sg + Sm 

^ JM ^ JM 

We use a system of units for the gravitational constant G and the light speed c 
such that = 1. The spacetime M is four-dimcnsional and is endowed with two 
conformally-related metrics: the Einstein metric denoted by g/j^^, and the Brans-Dicke 
(or physical) metric denoted by 

In the latter, a{ip) > is a coupling function entering the theory, and one recovers 
classical general relativity by choosing a{ip) to be constant. 

This theory admits two scalar fields: one of them, <p, represents the matter content 
of the spacetime and the other, tp, generates the gravitational field. When a is not 
just a constant, the matter field 4> does interact with the physical metric ^^j/, whereas 
the gravitational field equations for g^^, ip are formulated in terms of g^v, only. 

2.2. Choice of coordinates 

Throughout this paper we use a notation consistent with the one in Christodoulou [7] , 
in order to make easier the comparison between the two models. We consider a general 
spherically symmetric spacetime whose metric is expressed in Bondi coordinates [3] as 



g = -e^'^dv? - 2e''+^dudr + r'^dn'^, 
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in which the metric coefficients A, u depend on the coordinates u, r, only, and rfO^ 
represents the unit round metric on the 2-sphcrc. The relevant components of the 
Ricci tensor R^r and Rgg arc found to be (see for instance [3, 24]): 

_ 2 fdX dv^ 



r \dr dr 

\or or J 

with Rip^ = (sinf?)^ Ree- Observe that these formulas involve first-order derivatives of 
the metric coefficients, only. For completeness, we also determine the other two non- 
vanishing components of the Ricci tensor, which now involve second-order derivatives 
of the metric, i.e. 

R = e"^^"- 




, r du dudr dudr ) 





dr dr \dr ) dr"^ } dudr dudr ' 
Note in passing the following simple relation: 



2dX 

r du 



Alternatively, one may consider the Einstein tensor Ga^ = Rap — \Rga.i3 for the 
metric (/a/3) and compute its essential components 



Gun = 2— + 2e^-^ — 

r \ dr r du J 

_ e^"-^^ /e^^-l dX 

Gur — 1" 2 — 

r \ r dr 
2 fdiy dX 
dr dr 



r 



2.3. Evolution equations 

By varying the action of the theory with respect to both metrics, one gets two 
conformally-related stress-energy tensors T"'^ and T"'^ = a~^{^)T°'^ . Denoting by 
T = QapT"^ and T = QapT"^ their traces, one can check that (cf. [11] for details): 

Ga/3 = 2T„^ + 2 (daipdpi^ - ^ga^g'^'d^^jd^A (2.1) 



and 



Tap = da(t> dpcf) - ^ ga0 g'"^dp(t) d„(f>, 

Tap = a^{i')fap = a\iP){da<t>df}(f> - ^ga0g'"'dpct>d^ct>). 



Existence of naked singularities in Brans-Dicke theory 6 
Following Christodoulou [7], we use the future-directed null frame (n, /), defined 



by 



ou or 

Or 



By contraction of the Einstein equations (2.1) with n, I, we obtain the following system: 
dX dv f dip^ ^ 



Or ' dr 'ydrj ' ^'''^ 

dr dr r ' 



^ = r ( —— - e^^'"^ ( ^ 
du \drdu \ du 

Note that these equations involve only first-order derivatives of the metric and scalar 
fields. 

On the other hand, by defining the function 

the evolution equation for the scalar field ip reads (cf. again [11]) 

□gV = -^{i')T = a'(V')a(V').9"^ 9^30, 

where is the wave operator associated with the Einstein metric. In our gauge, this 
equation is equivalent to 

_ 2 + + e^''-^) + 2 ^ + ^t^^^ 

\ drdu r du J \ dr^ r dr dr dr 



Finally, the equation for the matter field (j) is obtained from the zero-divergence 
law for the stress-energy tensor: _ ^ 

VaT"'^ = 0, 

where V is the covariant derivative associated with the physical metric. In terms of 
the Einstein metric, the zero-divergence law for the stress-energy tensor reads 

where V denotes the covariant derivative for the metric g. This equation is equivalent 
to 
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a linear equation for which, in our gauge, becomes 



\ drdu r du J 
= -2a(V) (e 



\ dr"^ r dr dr dr 

(v-X) d'tp dcf) d<j) dip dcf) di]j 
dr dr dr du du dr 



2.4- Case of interest in this paper 

For dcfinitcncss and simplicity, wc study the case in which <7{ip) = cr is a constant, 
which corresponds to Brans-Dicke theory [6]. Integrating (2.5) we get 

a(V') = ao e^*, 

where flo is a dimcnsionless constant (independent of ijj). Assuming that it does not 
vanish, this constant can be eliminated by re-defining the matter field </>, and so we 
end up with 

a(V) = e'^^, 

which is the choice made in the rest of this paper. 

Our equations can be easily compared with those used in other investigations of 
the Brans-Dicke theory. For instance, Liebling and Choptuik [21] have (denoting their 
variables with subindex LC) 

(p = 'ipLC, <^^ = -^^LC, 2ct^ = (47r)ALc, 

which imply = e^^'^^ . (Special care must be taken with the fact that, in [21], units 
with G = c = 1 are used so that various factors like 47r arise in their equations.) 



3. Self-similar assumption and the reduced system 

3.1. Essential field equations 

We now impose continuous self-similarity on the solutions, that is, 

jl^sQi^u = "^Qixu (3.1) 

for some (conformal) homothetic Killing field denoted by S^^. To work with self-similar 

solutions, it is convenient to use adapted coordinates, in which the integral lines of 
are now coordinate lines. Every spherically symmetric and self-similar spacetime 
(except Minkowski) has a singularity at a point on the central world-line, and we shall 
use it to define the origin of time, so that u = represents the future null cone of 
the singularity. Moreover, since we are interested in the process of the formation of 
singularities, we (principally) work in the past region u < 0. 
We define Bondi's self-similar coordinates by 

x=—, T = -log{-u), (3.2) 
—u 

where the sign in r is chosen so that u and r increase simultaneously toward the 
future. The homothetic vector is now S = —dr, with integral lines x =const., and 
points away from the singularity — which is now located at r = -|-oo. Referring for 
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instance [15] for the general geometry of self-similar spacetimes, we note that in these 
coordinates the metric reads 

g = ((-e^'' + 2xe''+^)dT'^ - 2e''+^dTdx + x'^dfl'^) , (3.3) 

where v and A arc functions of (r, a;). The Lie derivative Cs is now simply —d/dr 
and, consequently, the symmetry condition (3.1) implies that all metric coefficients 
depend on x, only, so 

V = vix), A = X{x). 

This implies similar conditions on the scalar fields and, since arises in an 
undifferentiated form, the relevant condition is 

d 

—tjj = 0, implying V = V'(a;)- 

However, <p only enters the equations in differentiated form and hence the condition is 
d 

—9^0 = 0, implying ^ = x{x) + kr, 

where X = Xi^) is a function of x, only, and fc is a dimensionless real constant. 

Remark. An interesting variant of the above symmetry assumption was adopted in 
[17], where self-similar, complex-valued, scalar field solutions are constructed from the 
ansatz e'"'^^(a;), for some dimensionless real constants. 



Following Christodoulou [7], we define 

/3 = 1 



2x ' 

e=x^{x), (3.4) 
^ = xil)'{x), 

and we emphasize that (3 will replace v from now on. The Einstein equation (2.3) 
becomes 

and, by adding (2.2) and (2.3) together, we find 

2x^=e + e'"^^' - (e'^ - 1). (3.5) 
dx 

The third Einstein equation (in combination with the other two equations) yields the 
constraint equation 

e^^ = 1 + e e^"^ + {{6 + kfe^''^ + e) ■ (3.6) 

1 jJ 

Moreover, the wave equation for the matter field reads 
f}f} 

px— + (1 - (1 - p)e^^) e + k = -a{k + 2(59)^, (3.7) 
and the wave equation for the Brans-Dicke field is 

px$- + (1 - (1 - /3)e2^) ^ = ae^^'^ik + me. (3.8) 
dx 
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(We emphasize that there is no factor 2 in the parenthesis of the right-hand side of 
this last equation.) 

In conclusion, using the constraint (3.6) to eliminate A, we arrive at a system of 
four ordinary differential equations: 

d/3 



^ dx 



1 _ yt2g2a^ _ ^^2ai, ^ 2k)0 + + l) /3, (3.9) 

px^ = k{ke^'''''9 - 1) - a{k + 2/36»)^ + (e2'^'^(6i + 2k)9 + - l) p9, 

px^ = k^e^""^^ + ae^'^'^ik + pe)0 + [e^'^'^ie + 2k)e + - l) 

x^ =i 
dx 

The first two equations above reduce to Christodoulou's equations (see (0.27a) and 
(0.27b) therein) provided a and ^ vanish identically. Note that we can simultaneously 
change the signs of 6 and k without changing the structure of the system. Hence, 
without loss of generality, we can assume that A; > 0, while 9 still can have any sign. 

3.2. Reduced system 

We have found it convenient to rescale 9 with an exponential factor e"^"^, which 
compensates for the discrepancy by a factor 2 between equations (3.7) and (3.8) (as 
pointed out earlier). We define 

e = e^^e, 

and the remaining exponential terms can be combined with k into a single variable: 

K = e'^'^fc. 

The notation is intended to compare with Christodoulou's case, for which Q and K 
coincide with 9 and k, respectively. Using the variable s = Ina; and denoting d/ds as 
a prime, the system (3.9) now reads 

f3' =i-K^ ~ {2Ke + + + 1) 

/3 e' = K{Ke - 1) - cf{k + i3e)^ + {2kq + + _ -^^ /3e(3.io) 

Pi' = K^i + a{K + pe)e+ {2KQ + 62+^-1) /3^, 
K' =aKi. 

It is also convenient to make the change of variable 

1 

such that the system under consideration becomes polynomial in all variables: 
a' = {aK'^ + 2K@ + 6^ + + l) a - a^, 

e' = [aK^ + 2KQ + e'^+f-l)Q- a{aK + 9)^ - aK, (3.11) 
e' = {aK'^ + 2K@ + - l)i + a{aK + 6)9, 

K' = aK^. 

From now on, we refer to these equations as the reduced system, which is our main 
object of study. We sometimes use it in the form (3.10), evolving (3 instead of a. 
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Observe the combination 

aK'^ + 2KQ + Q^ + ^^ = {a- l){e^^ - 1). (3.12) 
We will later use the variable 

that obeys the evolution equation 

^{L^y = {aK^ + 2KQ + L'^ - l) L'^ - aKQ. (3.13) 

Remark. From O and ^ one can form a complex function A = Q + i^, with norm 
|A| = L, which satisfies the differential equation 

(7 \ K^7 

with Zq = K + ia and Z\ = 2K + ia. 

4. Interior solution originating at the center 

4-1- Analytical strategy 

The construction of our spacetimes, as solutions of the reduced system (3.11), will be 
performed in several steps; a main difficulty stems from the fact that the equations 
become singular for several values of the x-coordinate. This happens at the center of 
spherical symmetry, a; = 0, and at the self-similarity horizons corresponding to those 
values of x for which the homothetic vector S becomes null. In our Bondi coordinates, 
this corresponds to the condition g{dr-idr) = grr = for the past light-cone, which 
reads (see (3.3)) 

-e"^" + 2xe''+^ = 0, or /3 = 0. 

This section describes the construction of the past of the singiilarity, namely the region 
between the center worldline and the past lightcone of the singularity, the first self- 
similarity horizon. Following Christodoulou, we shall refer to this region as the interior 
solution. 



4-2. Critical points and regularity at the center 

We begin by determining all critical points corresponding to equilibria of the reduced 
system, provided cr ^ 0. In view of the right-hand side of (3.11) and provided 
the unknown functions a, Q, K have vanishing derivatives, only the following 
alternatives can arise: 

(a) a = 0, 9 = 0, ^ = 0, ii" arbitrary, 

(b) a = 1, e = -K, ^ = 0, K arbitrary, 

(c) a = 0,e = 0,^'^ = l,K = 0, 

(d) a = 2,Q = 0,^'^ = l,K = 0. 



Existence of naked singularities in 



Brans-Dicke theory 



11 




Figure 1. Schematic representation of the different regions in which the the 
solution is divided, matching at the past hght-cone (x = x*) and the future hght- 
cone {x = +00). The dashed lines denote Bondi coordinate lines. 



Points (c) actually belong to the exact solutions a = 0, 8 = 0, ^ = ±1, K = Kqc^'^^. 
Note that this collection of fixed points is not an extension of Christodoulou's result. 
This is due to the fact that the condition ^ = is not preserved by the evolution. In 
fact, several fixed points in Christodoulou's problem are no longer fixed in our case. 
This is clear in Figure 2, which shows a projection of the evolution flow on a slice 
(O,^) of phase space. 

The center of symmetry must correspond to one of the above cases. We 
need to rely on physically motivated regularity requirements, to select one of them. 
Specifically, we impose that the center of symmetry is regular before the formation of 
the singularity. Recall that the Hawking mass is determined from the metric coefficient 
A by the relation 

1 = e (4.1) 

r 

In order to avoid a singular behavior at the center, we impose that the mass tends to 
zero, which implies A = at the center. The equation (3.12) then selects the fixed 
points (a) and (b), above. 

Regularity at the center also requires that v(x = 0) is finite, to avoid a coordinate 
time singularity. However, the value of v at the center is gauge-dependent and, by 
normalizing // to vanish at the center, the coordinate time u in (3.2) coincides with the 
proper time of the central observer. Consequently, /? in (3.4) behaves like /3 ^ —l/(2x) 
and, equivalently a ~ — 2:r at the center. In particular, this condition implies that 
a = at the center, which is consistent with the critical point (a) above, only. 

In summary, we obtain the following values for the critical point of interest: 

a = 0, 6 = 0, K = Kq, ^ = at the center a; = 0, (4.2) 
where Kq is an arbitrary non- negative constant. 
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Figure 2. Projection onto the plane {a = 0, 0, ^, K = 1/2) of streamlines of tlie 
flow vector field of system (3.11) for a = 1/3. We have marked the fixed point 
(a) at = 0,5 = and the projections at © = 0, ^ = ±1 of the exact solutions 
mentioned in the text. We have also marked the points = —K it V K'^ + 1, 
which are fixed points in Christodoulou's system, but are no longer fixed points 
in our system. Those four points now form part of an unstable projected structure 
that resembles a limit circle, though we do not know whether a true limit cycle 
exists in the full phase space. 



4-3. Linear stability of the critical point at the center 



After linearizing around the critical point (a), the Jacobian matrix of the system in 
the linearized variables {5a,6Q,5^,SK) reads 



/ 1 ^ 

-Ka -10 

0-10 

V aKo / 

Its eigenvalues are 1,-1,-1,0, with respective eigenvectors 

,(1) - -Kof^ (2) _ 1 

' 



V 
/ 

1 



.(4) ^ 



V / 




V 1 J 



respectively. (Figure 2 corresponds to the two negative eigenvalues.) The origin (a), 
therefore, has an unstable branch, tangent to the vector e^^\ and the corresponding 
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e 

V K J 



( \ 





V ) 



+ ai p/ 



V 



1 






+ o{e% 



where a\ is a parameter. 

Imposing tlie normalization ax~^ —2 at (a), we get a\ = —2 and we obtain 
an interior solution in the neighborhood of the center, satisfying 







\ 


1 


-2e" 


\ 




e 























v 


K 


I 


V 




I 



(4.3) 



We now show that this imphcs that the spacetime is (mildly) singular at the 
center. We start by taking the trace of (2.1), which gives the Ricci scalar 

i? = _2r + 2(5"^a„^a^^). 

Taking the form (3.3) of the metric, wc get 



R = 2e^ 



-2e 



{a'^dr<t>dx<t> + drtpdxtlj] 



(4.4) 
(4.5) 



+ e 



-2A 



2xe 



(9,0) V (5, V)')) 



Using the reduced system (3.11) together with the expansion (4.3), we get (for the 
first two variables) 

a = -2e'- 46^^ + o (e^^) , 
e = i^oe^ + ^e2« + o(e2«). 
Therefore, the local expression for the Ricci scalar becomes (we have replaced e* with 

X) 

R = e'^^ {-2K^-4K^x + o{x'^)), (4.7) 

where we have used limx^o i^i^) = hm^j^o M^) = 0- 

Consequently, the presence of a non-vanishing linear term —AKqX in (4.7) shows 
that, when viewed as a geometric object in a spherically symmetric spacetime, the 
scalar curvature R is continuous but not difFerentiable at the center x = 0. (Only 
even powers of x should, otherwise, be allowed.) Hence, the spacetime contains a 
(mild) singularity before the central curvature singularity forms at u = 0. 



4- 4- Integration in the interior region 

We denote by (— oo, s*) the maximal interval on which the solution is defined, and 
now check that, provided s* is finite, a must blow-up. 

Claim 1. Either the solution is defined and regular up to s* = +cxd, or else 
a{s) — )• —DO but 0,^ and K remain hounded as s ^ St,. 
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Indeed, sc suppose that is finite and we are going to prove that both Q and 
^ are bounded. First, a < for all s < s*. Indeed, a{s) < at least for an interval 
of the form (— oo, sq] with sq < «*• Thus, /9(s) = l/a(s) is negative and finite on the 
same interval. Also, setting B — 2KQ + 0^ + + 1 and C = 1 — K'^, from the first 
equation of (3.10) we deduce that (for so < s < s*) 

I3{s) = /3{so)e- -^4 + r C{u)e- 

J so 

Thus, I3{s) remains finite on (— oo, s*) so that a{s) = 1/I3{s) can not vanish in this 
interval and remains negative. Now, the equation (3.5) reads 

2A' = + ©') - (e'^ - 1) 

thus we have 

r_^e^'9{s')ds'' '^'^ 

But, L can not vanish, except at an isolated point. Indeed, if there is si < s* 
such that L{si) = 0, then 6(si) = with Q'{si) = —a(si)K{si) > 0. Thus, functions 
© and L are both different from zero in the neighborhood of si . We obtain that the 
function g is strictly monotone increasing, 

r e''g{s') ds' < g{s) f e«' ds' = e'g{s), 

J — oo J — oo 

so that e^^ > 1, i.e., A(.s) > for all s < s*. Setting 

T = aK^ + 2eK + L"^, 
and using (3.6), we obtain 

e^^ - 1 = 

a — 1 

Since a < for s < .s*. wc get r(s) < for all s £ (— oo, s*). 
Let us introduce now the quantity 

H = +f + 2- = + 2-. 

a a 

Using (3.13) and the third equation in (3.11), we obtain 

H' = TH + TL^ - 2(^2 + -). 

a 

Then, completing the square in the last term, 
H'<TH^^. 

Now, fix So < s*. Using the Gronwall inequality, we get for s G [sq, s*), 
H{s) < H{so)e^'o ^(^^"^ + ^ f e^" ^^^^''^du 

J so 

< H+{so) + ^(s - So) < H+{so) + ^(s* - So), 



Existence of naked singularities in Brans-Dicke theory 15 

where H~^{so) = max(0, i?(so))- We deduce that both and ^ are bounded, and 
since K{s) = K{so)e-'''o^ , K is also bounded. We conclude that necessarily a — oo, 

cLS S ^ S ^ • 

Claim 2. Assume that s* < +oo. Then, there exists a real Q < < 1 such that, 
lim K{s) = K„ 

lim e(s) 



a 



lim ^(s) = - „2 , 2 • 

To estabhsh this claim, we proceed as follows. First, since ^ remains bounded, 
I-oo ^(^)^^ = where S is some real constant, so that 

hm K{s) = Koe"^ = K^, 

S^S* 

with > 0. Also, the first equation of (3.10) gives ^ 1 — K^^ and we obtain 

if*' - 1. (4.8) 



But since /3{s) < for ,s < ,s», we get that necessarily Kf, < 1. The case iV'* = 1 will 
be excluded at the end of the argument. 
Now, let us introduce the new variables 

5i=Ke-a^- 1, S2 = K^ + aO. (4.9) 

Thus, the expressions of ©' and ^' in (3.11) become 

e' = j6i + hi, ^' = j62 + h2, (4.10) 

where 

hi = {2KQ + + C' - 1)9 - (7^9, 
/i2 = + + ^2 - 1)^ + ctO^ 

Using (4.9), (4.10) and the last equation in (3.11), we get 

5[ = j{K5i-a52)+h, 5'2 = j{(j5i+K52)+l2, (4.11) 

where 

li=(jKiQ + Khi-ah2, (4.12) 

h=(TKS^ + Kh2+ahi. (4.13) 

Thanks to (4.11), the quantity 5 ^ {5l + (5^)^/^ satisfies 

6' = ^S + l, (4.14) 

where 

hSi + 1262 
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By assumption, li and I2 are bounded on (— oo,s,) so that I is bounded too. 
Choosing s < s*, and integrating (4.14) over [sq,s] we obtain 

6{s) = e-^(") U{so) + f e«("') l{s')ds'\ , 

J Sq 

where 

Function ( is increasing, and by (4.8), it tends to +00 as s ^ s*, together with its 
derivative Thus, 

e-C(s) f e«(^')ds' as s s„ 

J So 

and since I is bounded we obtain 



Hm 5{s) = 0. 
Finally, we can rewrite (4.9) in the form 

and we get 

hm e(s) = „ , hm^(s) = -- 



Now, it remains to prove that if* < 1, i.e., if* = 1 is excluded when < +00, 
as weU as if* > 1. So, assume by contradiction that if* = 1. Thus, the last equation 
in (3.11) reads 

2 

K{s) = 1 - (.s - 5*) + oi\s - s,|), (4.15) 

Also, (3.10) gives /3'(s) — >■ as s — >■ s*, so that /3(s) = o(|s* — s|). More precisely, 
using (4.15) in the first equation of (3.10) we obtain the following expansion 



I3'{s) = 2^—{s - s*) + o{\s - s*|). 
a + 1 



Thus, P is decreasing in the neighborhood of s* which is impossible since /3 < and 
/3(s) — )• as s ^ s*. 



4-5. Conclusions for the interior region 

In the interior region, the situation is similar to that of Christodoulou, in the sense 
that the presence of a first self-similarity horizon, the past light-cone of the singularity, 
is determined by the value of a constant if*, which must be below the threshold 1. The 
main difference is that in Christodoulou's case this constant is the parameter k, fixed 
throughout the problem, while in our case the constant if* is dynamically determined 
by the evolution and, so, depends on the set of initial conditions, especially the initial 
value ifo of the variable if. 
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Figure 3. Evolution of K(x) starting from values Kq ranging from to 7 in 
steps 1/2, hence 15 curves. The coupling constant is cr = 3/2 in all cases. Each 
evolution stops when reaching the singular point a;,, which for Kq = corresponds 
to Minkowski, with = 1/2. For reference, a dashed line represents the value 
K = 1. 



We have performed numerical integrations of the reduced system of equations 
to investigate how K evolves. The results are summarized in Figure 3, which shows 
several evolutions of K starting from different values at the center. In all cases we see 
a decrease of K until values below 1 are reached, and then we reach the singular point 
X* = e"' , where the integration is stopped. We have found this behavior in all tested 
cases, including cases with large values of the initial Kq (above 1000, say). The value 
of a does not alter the qualitative picture, though the decay of K is faster for larger 
values of a. Interestingly enough, for large values of Kq the final value K^, is almost 
independent of that initial value. For sufficiently large values of a we find numerically 
that tends to 4/(3(t). 

In other words, we can have solutions in which the central singularity has a past 
light-cone for initial values of the constant k for which Christodoulou's corresponding 
solution would be more pathological, with that light-cone becoming actually a border 
of the spacetime. 



5. Extension to the exterior region 



5.1. The singular points 

According to the previous section when s, < +oo, the solution of (3.11) converges to 
a singular point of the form {(3,Q,£_,K) — (0, 0*, A',) — Uk,, where < A'* < 1, 
0* = K^-^a'^ ' ^* ~ ~~ A'^+o-^ • treat the solutions in the neighborhood of such 
singular point, and following [7], we introduce a new independent variable t satisfying 

dt ^' 

which converts the singular point at finite s into a critical point at t — > -co. 
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Thus, by using (3.10), the variables /?, O,^, and K satisfy the system 

dp 
Itt 

^ =pQ{2K@ + Q'^+^''-l-aCi+K'^Q-K-aK^, (5.1) 



^ = (1 - K^)p - (2X9 + + f + l)/32, 



dt 



dK 



The singular point Uk, is an equilibrium point of the previous system, and the 
Jacobian matrix at this point reads 

( l-Kl \ 



Ki{Ki+i) 

(/<-J+a2)2 



Kl -aK, 



\ 



The spectrum of A{Uk,) is given by 
Sp{A{Uk^) = {Q, 1- 



iuK^, Kl + iaK^}. 



The eigenvalue corresponds to the fact that the set C* = {Uk,,0 < -f^* < 1} 
defined by the equilibrium points of the form Uk,, with < if* < 1, is a (one- 
dimensional) curve. Each point Uk, has an unstable manifold W^* of dimension three, 
corresponding to the three other eigenvalues having a positive real part. Naturally, 
Wf^' must be transverse to C*. 

We also observe that all solutions originating at Uk, and extending from the 
interior to the exterior region, admit the following expansion (when t — )• — oo): 

/ ^ \ 

e 



\K J 



Uk, + aie(i-^-)*Ci + a2e^-*C2(0 + a^e^^'^Cait) 



+ o(e 



Kit p{i-Kl)t 



), 



where C2{t) and C3{t) are (bounded) periodic-rotating vector- valued functions, and Ci 
a fixed eigenvector of the matrix A{Uk, ) corresponding to the eigenvalue Ai = 1 — K^. 
The three vectors Ci, C2{t) and C:i{t) arc linearly independent for all t. Up to a 
translation in t, we can assume that ai = 1 and, thus, we obtain a two-parameter 
family of solutions. 

The expansion above shows that the functions on the left-hand-side are continuous 
at the past light-cone, but not differentiable in s there for < < 1, since they will 
generically contain terms of the form (s — s*)^* / (i-^.) or (s — s^.)^^^^^' . In particular 
will be only continuous, and hence curvature will be discontinuous, though still 
finite, at the past light-cone. 
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5.2. Exterior solutions 



We computed in Section 4.2 the fixed points of our dynamical system. In the interior 

the relevant point was (a), but now we will study (b) and (d). We first focus our study 
of the exterior region on an equilibrium point of the form (b), namely 

Pk. = {1,-K*,0,K*) for a fixed < K* < 1. 

The Jacobian matrix of system (3.11) at this point is given by 



The spectrum of this matrix is given by 

Sp{A{Pk-)) = {0, K*^-l, -1-iaK*, -1 + iaK*}, 
which gives the asymptotic behavior 



/ 


K*"^-! 








\ 




-K*^-K* -1 





-1 






-crK*'^ -aK* 


-1 


-aK* 




V 





(jK* 





/ 



/ a 

e 



= Pk* + aie^'''"-^>Ci + a2e(-i-'^^*)*C2(t) 
+ a3e(-i+'^^*>C73(s)+o(e- 



(5.2) 

The eigenvalue corresponds to the fact that the set C* = {Pk*,0 < K* < 1} 
defined by the equilibrium points of the form Pr-. , is a (one-dimensional) curve. Each 
point Pk* has a stable manifold of dimension three, corresponding to the three 
other eigenvalues having a negative real part. Naturally, must be transverse to 
C*. 

Consider now the two isolated critical points of form (d) 

K = (2,0,e,0), e = ±l 
The Jacobian matrix of system (3.11) at Vc reads 



( -2 



V 



4e 

-o-e 

2 





\ 
-2(ae+l) 



ae ) 



The spectrum of the previous matrix is given by 

Sv{A{V,)) = {-2, -a, <7, 2}. 

Each point V^, e = ±1 has a stable manifold w'f^ and an unstable manifold w/^\ 
both of dimension two. 

The points Pk* are attractors (except for the marginal direction connecting 
them), and our numerical simulations below will show that it is indeed possible to 
evolve from points Uk, at the past light-cone to points P^* at the future light-cone. 

In this section, in order to establish that this is indeed the future light-cone 
u = of the singularity, we investigate the behavior of exterior incoming null rays for 
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solutions terminating at such point Pr-. . The exterior condition means that we work 
with 13 > 0. Evolving towards a fixed point means that wc can approach x = +oo, 
and hence this is either r = +oo or u = 0. Therefore we need to show that incoming 
null rays can reach x = +oo at finite r. In self-similar coordinates (r, s) the equation 
of incoming null rays is 

- = -V\ - = /3, 

du 2 ' dr ' 

and hence, for a given ray originating at (tq, sq), 

T — To = / a{s')ds' . 

J So 

Prom (5.2) we get that 

{a{s) - l)e(^~-^*^'* c as s +oo 
for some constant c. Therefore, 



log 



r = s — T = so — Tq — / {a{s') — l)ds' 

J sn 



converges to a finite quantity as s — )■ +oo. 

Finally, it remains to show that the future light-cone is not a curvature singularity, 
so that the spacetimc can he continued beyond it. Indeed, while the self-similar 
coordinate system (r, x) becomes singular on the future light-cone, we still can take 
the limit x +oo in the formula (4.6), replacing the results in (5.2). The result for 
the Ricci scalar R of the spacetime metric is 

r 7T— ^ — as a; — >■ -|-oo. 

{1 + K*^){K*^ + a^) 

Hence R is Unite everywhere on the future light-cone — except of course at r = 
(which is a curvature singularity). Similar expressions can be derived for the Gauss 
curvature of the two-dimensional reduced spacetime, or for the Kretchmann scalars of 
the four-dimensional and two-dimensional metrics. 

Our points Pk* = {l,—K*,Q,K*) play the same role as Christodoulou's point 
Pq = (1, —k), and we see that they are indeed closely related since the field ^ vanishes 
at our points Pk* ■ This can be interpreted as a sign that the Brans-Dicke field is 
becoming irrelevant on the future light-cone of the singularity and, therefore, that the 
spacetime has the same properties as the ones of Christodoulou's solutions within a 
small neighborhood of the light-cone. 

An important difference between our construction and Christodoulou's one is 
that he can construct time-symmetric solutions, since the data on the past and future 
light-cones coincide. But, such construction is not possible here precisely because the 
point Uk, has a non-vanishing field ^ while the point Pk* has ^ = 0. 

Christodoulou uses the above fact to "copy" the past region of the spacetime onto 
the future region (cf. Fig. 1), finding a possible complete spacetime containing the 
naked singularity, with a center which is only mildly singular. In our case, to have a 
complete spacetime, we would need to evolve (numerically) further from data on the 
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future light cone — but we shall not do that here, as we decided focus on establishing 
the presence of the singularity, only. 

Finally, the structure of phase space around this point is quite different in our 
system. There is no equivalent of Christodoulou's points P+ and P_, due to the 
general tendency of the variables Q and ^ to "rotate" among them (which is the 
same phenomenon as the one discussed in Figure 2). The two-dimensional funnel in 
Christodoulou's pictures is converted here into a higher-dimensional analogue. In the 
following section we use numerical evolutions to demonstrate this behavior. 

6. Numerical investigations 

6.1. Numerical strategy 

We perform a numerical integration of the system (3.11) using a fourth-order Runge- 
Kutta scheme. This is done in four steps: 

(i) Using the coordinate x = e^, from x = 0, with a = 0, Q = 0, K = Kq, and ^ = 0, 

we integrate the equation (4.2) up to some finite value xq = e^". 

(ii) Then, switching to the reduced system (3.11) in the coordinate s we use an 
adaptive-step approach to integrate on the interval [so,s*) and get as close as 
possible to s*, where lims-,.^^ a = — oo. 

(iii) Next, starting again from si > s^, wc integrate the system (3.11) backward to 
reach s*. At Si, we have two new parameters, namely Qi = 6(si) and = £,{si), 
whereas a{si) and K{si) are determined so that all quantities (^, O, ^, K) match 
at s = 5 ^ . 

(iv) Finally, from s = si, we integrate forward with a constant step-size and check 
whether the system (3.11) converges to a stationary point or diverges. 

A first step is needed to start from the accurate values at the central singularity, as 
defined in Section 4.2; the matching with the second step at s = Sq is straightforward, 
since the transition is done at a point at which all quantities have regular behavior 
and we only make a change of coordinate from x to s. The matching at s* is much 
more complicated and we first explain now the technique used to recover the results 
by Christodoulou [7]. 

6.2. Numerical integration of Christodoulou's system 

In his study, Christodoulou solves the equivalent of our two first equations (for 
a = a{s) and 6 = 6(s)) in the reduced system (3.11): 



(Cf. the equations (1.1a) and (1.16) in [7].) In a neighborhood of s,, for s > s*, 
Christodoulou finds that the solution 6 depends on a real parameter ai (see (2.13) 
and (2.14) in [7]), but one always has 




(6.1) 



(6.2) 



lim 6 = —. 
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!2 



k=0.3 
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= 0.5 
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= 0.5 




s 



Figure 4. Solutions = a{s)~^ and 9{s) in the case of Christodoulou's 

system (6.1)-(6.2). The matching of the solution at s = s* is shown for two 
different values of the parameter 9i = 9(si). 



On the other hand, the solution /3 — a ^ has the foUowing behavior: 

/3 = (1 - k^) (s -s,) + (\s ~ s,|') . (6.3) 

With these information, we devise the numerical integration strategy as follows. Given 
a value of A:, we integrate the system (6.1)-(6.2) until s — > s*, and thus determine the 
approximate value of up to high accuracy. We then define 

si = (l±e)s*, (6.4) 

with the sign chosen so as to s\ > s, and e ^ 0.03 for numerical convenience. We 
then set 

"^^^^= (1-P)e.. ' 

e{si)^9u (6.6) 

with 9i a new parameter than can be freely chosen and is to represent the degree 
of freedom, induced by the parameter oi of Christodoulou's study (see Section 2 
and Eqs. (2.13b)-(2.14c) of [7]). From the two points (si, ai), (si, ^i), we integrate 
backward toward s* and a{s) — ^ +oo. When doing so, we find numerically that, in 
most cases, a is diverging at some value s'^ 7^ s*. This is due to the approximate value 
of Qfi in (6.5) that was chosen to initiate the integration. Since we are dealing with 
an autonomous system, we can perform a slight shift As — s'^ ~ s* in the variable s, 
so that lim^^^+ a{s) = +00. 

In Figure 4 are shown the numerical solutions of the differential system (6.1- 
(6.2), for Kq = 0.3 and with two different values of 9i = 0{si). We have numerically 
observed that, if the parameter e was small enough and for any value of 6*1, the 
backward integration detailed here-above would always bring back to the solution 




^(s*) = and 9{s^) — 1/k. We were thus able to numerically recover the result 
by Christodoulou [7] that, for a given A: < 1, each solution of the system (6.1)-(6.2) 
connects to a one-parameter family of solutions at s = s*. 

By varying the parameter 6i, one can reach different asymptotic regimes, when 
s —7- +oo: both fields a and 9 can diverge or they can converge to the stationary point 
Vq : {a — 1,9 = —k). This last case is displayed in Figure 5, with the trajectories of 
the solutions in the {9, a) plane. (This is to be compared with Figure 4 in [7].) 

6.3. Numerical solutions of the reduced system 

In the case of interest in this paper, we have to deal with the matching of four fields 
{a,Q,^, K) at s — s^. We use the same numerical technique, with the four steps 
described at the beginning of this section. We have numerically observed that here, in 
addition to 0i(6i = 6(si)), we need to specify the value of = $(si). On the other 
hand, near an equilibrium point Uk,, with < A'* < 1 we get from the first and last 
equations in (3.10): 

so that we can write the following expansions: 

/3(s) = {l-K^,){s-s,) + o{\s-s,\), (6.7) 

K{s)^K,.- ^\ (s- s,) + o{\s- s,\). (6.8) 

We can therefore obtain numerical estimations of the values of these two fields at 
s = si and perform the integration backward, from si toward s*. We again do the 
shift in s to match f3 up to machine precision at s = s*, but doing so does not allow 
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Figure 6. Solutions /3 (= a ^) ,0,^ and K of the reduced system (3.11). The 
matching of the fields at s = s« are shown for three sets of the parameters 

(Gi = 0(si),ei = e(si)). 



for an accurate matching of K{s). We then do several (usually no more than five) 
integrations from si toward s*, correcting each time the starting value K{si) in such 
a way that, at the end, the function K{s) is continuous at s*, up to machine precision. 

Results are displayed in Figure 6, where we have taken Kq = 0.3 and a = 0.1. The 
matching of all four fields are done with the setting of two new parameters (0i,^i), 
which seems to indicate that a solution starting from the interior region (s < s,) 
connects to a two-parameter family of solutions in the exterior region (s > s,). 
Depending on these two parameters (Oi, ^i), we numerically recover a behavior similar 
to that of Christodoulou's system [7]: for some values of k, cr, 9i and ^i, the system 
can converge to the stationary point (b) of Section 4.2, that is, (1, —K*, 0, K*) with 

K* = lim K{s). (6.9) 

Part of this behavior is displayed in Figure 7, where the trajectories in the {9, a, ^) 
space, for three sets of parameters (9i, fi). Each set of parameters can lead a priori to 
a different limit Pk* ■ If these trajectories were projected onto the (a, 6) plane, they 
would resemble a lot the ones of the general-relativistic system of Figure 5, studied by 
Christodoulou [7] , with the noticeable difference that we no longer have a single limit 
Vo, but the endpoint depends in general on the value K* (see (6.9)), which changes 
from one set of parameters to another. 

7. Conclusions 

We have studied the formation of naked singularities in the process of the gravitational 
collapse of a real massless scalar field, and have generalized Christodoulou's 
construction of a family of spacetimes containing a naked singularity. While 
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Figure 7. Trajectories of the solutions of the reduced system (3.11) in tlie {6, a, ^) 
space, for three different sets of parameters (0i,^i). 

Christodoulou worked within the classical Einstein theory, we have here considered 

the Brans-Dicke theory. This effectively led us to deal with a new scalar field and to 
promote Christodoulou's constant parameter fc to a dynamical variable K. 

We were able to fully analyze the interior region and rigorously justify the 
matching across the past light-cone of the naked singularity. Partial analytical 
information was also obtained in the exterior region, and we finally completed our 
study with numerical simulations. We could show that the variable K always decreases 
from any initial value to a value fC* smaller than 1 at the past light-cone. This 
eliminates the possibility in Christodoulou's work of forming pathological solutions 
without a past light-cone. In that sense the addition of the Brans-Dicke field has a 
"regularizing" effect. The behavior of the solutions in the exterior region is similar 
to the one of Christodoulou's solution and, in fact, the Brans-Dicke field vanishes on 
the future light-cone of the singularity (a Cauchy horizon). However, the dynamical 
structure of the phase space is quite different, as our phase space is much larger and 
does not contain Christodoulou's case as a subspace. 

Like Christodoulou established for solutions to the classical Einstein equations, 
our solutions are probably highly unstable against small (radially symmetric, for 
instance) perturbations. 

It is quite reasonable to expect that more general scalar-tensor theories would 
exhibit a similar behavior. In particular, it would be interesting to extend our 
conclusions to the more general models arising in the so-called f{R) theories of gravity 
when the action involves a nonlinear function of the Ricci scalar. 
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